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Abstract 

We describe fluctuating two-dimensional metallic antiferromagnets by transforming to a rotating 
reference frame in which the electron spin polarization is measured by its projections along the local 
antiferromagnetic order. This leads to a gauge-theoretic description of an 'algebraic charge liquid' 
involving spinless fermions and a spin 5 = 1/2 complex scalar. We propose a phenomenological 
effective lattice Hamiltonian which describes the binding of these particles into gauge-neutral, 
electron-like excitations, and describe its implications for the electron spectral function across the 
entire Brillouin zone. We discuss connections of our results to photoemission experiments in the 
pseudogap regime of the cuprate superconductors. 



I. INTRODUCTION 



An understanding of the nature of the electron spectral function in the underdoped 
'pseudogap' regime has emerged as one of the central problems in the study of the cuprate 
superconductors. A wealth of data has appeared in photoemission studies, some of which™ 
has been interpreted using a model of "Fermi arcs" across the Brillouin zone diagonals; 
other studies^! have indicated the presence of pocket Fermi surfaces in the same region of 
the Brillouin zone. Scanning tunnelling microscopy (STM) studies^ also indicate a Fermi 
arc of excitations which appears to end abruptly at the magnetic Brillouin zone boundary 
associated with two sublattice Neel order. Another issue of interest in experiments has been 
the angular dependence of the electronic excitation gap energy in the superconducting state. 
A 'dichotomy' has been notecP^^ between the behavior of the gap near the nodal points 
on the Brillouin zone diagonals, and the antinodal points along the principle axes of the 
Brillouin zone. 

Related studies have also been made of the electronic spectra in the electron-doped 
cuprates^OS. Here the spectra seem to be closer to that expected from conventional spin 
density wave theory, and so constitute an important testing ground of theoretical ideas. 

The recent observation of quantum oscillations in a strong magnetic fielcP^l has given 
further impetus to the development of a theory the normal state of the underdoped cuprates. 
Some of the transport dat d 23 * 25 ^ has been interpreted in terms of the presence of electron 
pockets in the hole-doped cuprates, for which there is no apparent evidence in the photoe- 
mission studies; the latter, however, have only been carried out without a magnetic field. 

Finally, an important motivation for our study comes from numerical studies^^ of 
Hubbard-like models for the cuprate superconductors. These studies show a significant 
regime without any antiferromagnetic order, but with a dichotomy in the normal state 
electronic spectra between the nodal and anti-nodal regimes. It would clearly be useful to 
have analytic effective models which can be used to interpret the numerical data, and we 
shall propose such models here. 

The theory we present here builds upon the framework set up in Refs. |33]and|3H These 
papers described a non-Fermi liquid state which was labeled an 'algebraic charge liquid' 
(ACL); this state was obtained after the loss of antiferromagnetic order in a metal via an 
unconventional transition^ . As we will review below, in its simplest realization, the degrees 
of freedom of the ACL are spinless fermions and a S — 1/2 complex boson z a which interact 
via an emergent U(l) gauge force. This gauge force has strong effects even in possible 
deconfined phases, and can lead to the formation of electron-like bound states between the 
fermions and z a . The main purpose of this paper is to present a general discussion of the 
dynamics of these fermionic bound states, and their influence on the photoemission spectra. 
A full description of the pseudogap regime in the hole-doped cuprates will also require®' 
considerations of pairing of these fermionic bound states into bosonic Cooper pairs: this we 
defer to a subsequent paper. 
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In our previous workP^l we presented analytic arguments based upon the structure of a 
continuum theory valid at long wavelengths. Here, we shall present a more general formu- 
lation of the theory, which allows computation of the electronic spectrum across the entire 
Brillouin zone. We shall show how arguments based upon symmetry and gauge invariance 
allow construction of an effective theory for the electronic spectrum. The theory will contain 
a number of coupling constants, whose values will have to be determined by comparing to 
numerical studies or experiments. Also, while the previous workP^l U sed strong-coupling 
perspective, starting from the Schwinger boson theory of the antiferromagnet. It is possible 
to derive the results presented below also from this strong-coupling approach. However, we 
will choose to present our results by departing^ from the "spin-fermion" model, which was 
originally developed from the weak coupling expansion. 

Let us begin by defining the Lagrangian, C s f, of the spin-fermion mo del We consider 
fermions Cj a (a, f3 =t, 4) hopping on the sites of a square lattice. These are coupled to the 
fluctuations of the unit vector field (a = x, y, z) representing the local orientation of 
the collinear antiferromagnetic Neel order. We will restrict our attention here to antiferro- 
magnetic order at wavevector K = (vr,7r), although generalizations to other K are possi- 
ble. Throughout, we will freely make a gradient expansion of over spatial co-ordinates 
r = (xi,yi), focusing on the long wavelength fluctuations of the order parameter. However, 
the fermion fields Cj iQ , have important low energy modes at many locations in the Brillouin 
zone, and so we will not make any gradient expansion on the fermion operators. We have 
the imaginary time (r) Lagrangian 

C s f = C c + C\ + C n 

l~-c — ^ ] C icSPr ~ H) c ia ~~ ^ ] tij (^ c ia C ja + C joAa^ 
i i<j 

i 

Cn = hS £r ^ drna)2 + v2{drna)2 } ■ (L1) 

Here are arbitrary hopping matrix elements describing the "large" Fermi surface, fi is the 
chemical potential, A is the spin-fermion coupling, g controls the strength of the antifero- 
magnetic fluctuations, a a are the Pauli matrices, and the n a field obeys the local constraint 
^2 a (n a ) 2 = 1. Almost all previous studies of the spin-fermion modeP^have involved a per- 
turbative treatment in powers of the coupling A, along with resummations of this expansion. 
Here, we will not expand in powers of A, treating it is a coupling of order unity: instead our 
analysis are motivated by expansions either in the number of field components, or by small 
9- 

For sufficiently small values of the coupling g, the model C s f clearly has an antiferromag- 
netically ordered spin density wave (SDW) ground state with (n a ) ^ 0. We are interested 
here in the mechanism by which this order is lost as g is increased, and a metallic state with 
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no broken symmetries is obtained. In a recent paper 371 with others, we argued that there 
were 2 generic possibilities: 

(i) In the first case, there was a direct transition at a single critical g = g c to a Fermi liquid 
metal with a large Fermi surface. This transition has been examined in previous workP^, 
and is directly expressed in terms of fluctuations of the 0(3) order parameter n a ; the T > 
crossovers above g = g c have also been described 23 * 3 -^. 

(ii) The second possibility involved intermediate non- Fermi liquid phases before the large 
Fermi surface metal was reached at sufficiently large g. In this case, the 0(3) order parameter 
was parameterized in terms of the spinor z a by 

< = (1-2) 

The spinor field z a is the natural variable to describe the loss of magnetic order at g — g c , 
and the non-Fermi liquid phases above g c , and replaces the 0(3) order parameter n a . 

Our focus in the present paper will be on the second possibility. Part of our motivation 
comes from transport measurements^ 4 ^, which show an extended regime of non- Fermi liquid 
behavior as T — > in high magnetic fields. We shall describe the photoemission spectra 
at non-zero temperatures on both sides of g c . Our g > g c results are candidates for the 
pseudogap regime of the cuprates, and relate especially to recent experimental results of 
Meng et alW. 

To complement the approach taken in Ref. [33 here we will motivate our choice of the 
non-magnetic non-Fermi liquid phases by extending the theory^ of the loss of Neel order 
in the insulator at half filling. Although motivated in a theory of metals, C s f also contains 
(in principle) a complete description of the insulating states at half-filling. Crucial to the 
description of insulators, 1 41 ' 42 1 are 'hedgehog' point tunneling events ('instantons') in which 
n a points radially outwards/inwards from a spacetime point. These hedgehogs carry Berry 
phases: in C s f the Berry phases are expected to appear from the determinant of the gapped 
fermions integrated out in a hedgehog field for n a . All previous treatments of the spin- 
fermion model have neglected the hedgehog Berry phases, and this may well be appropriate 
under suitable circumstances in certain superconducting states.^ 

In this paper, we wish to focus on regimes and phases where the hedgehog tunnelling 
events are suppressed. In the insulator, the hedgehog suppression is a consequence of quan- 
tum interference from the hedgehog Berry phases, and leads to interesting new 'deconfined' 
phases and critical points.^ 44 ' 45 1 Hedgehog suppression is also possible in certain exotic 
metallic states (to be described below), where fermionic excitations near a Fermi surface 
lead to a divergence in the hedgehog action. ^ ' 34 ' 43 ' 46 ' 47 -) 

It is useful to begin our analysis by adapting the phase diagram of Ref. [33] describing the 
doping of an insulating deconfined critical point - see Fig. [TJ In the insulator, at x = 0, the 
transition is between a Neel state and a valence bond solidP^This transition is described^ by 
a deconfined critical point, with the order parameter n a replaced by a complex 'relativistic' 
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Confining metal 
with possible 
broken square 
lattice symmetry 




FIG. 1. Proposed phase diagram^ of a quantum antiferromagnet doped with carrier density x. 
The strength of antiferromagnetic fluctuations in the insulator is controlled by the coupling g. At 
x = 0, there is a deconfined critical poinpSl, separating a Neel state and a valence bond solicP^. The 
broken symmetries in these two states survive in metallic Fermi liquid state at non-zero x. More 
importantly, the deconfined critical point broadens into a non-Fermi liquid phase, the algebraic 
charge liquid (ACL), which is the focus of attention in the present paper. The conventional large 
Fermi surface Fermi liquid is not shown in the phase diagram above. It appears at larger x, and 
its phase transitions to the ACL and the SDW states were discussed in earlier workPH 



boson z Q which carries charges under an emergent U(l) gauge field A^. this mapping will be 
reviewed in Section [IT} Upon moving away from the insulator, the fermions c a are replaced 
by fermions ip p fermions which do not carry spin, but do carry charge p = ±1 of the U(l) 
gauge field A p . We will present the Lagrangian of z a , ip p , A p which described the phase 
diagram in Fig. [T] in Section [TTJ 

In the SDW phase of Fig. [TJ the z a , if) p , A p theory is formally identical to the spin-fermion 
model of c a , n a . This is because hedgehogs are strongly suppressed when there is magnetic 
order. We explore this connection in Appendix [XJ and in Section IV to obtain the spectral 
functions of the electrons in the 'renormalized classical' (RC) regime®^. 

The transition out of the SDW phase in Fig. [T] is into a fractionalized metal phase, which 
was dubbecP^ an algebraic charge liquid (ACL). This transition is in the 0(4) universality 
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classic To leading order, the elementary excitations of the ACL are simply the constituent 
fractionalized fields z a , ijj p} A^, and their interactions have been discussed in previous workP^ 
Here, our primary results will be on the fermionic spectrum over a large intermediate 
length scale near the boundary between the ACL and the confining metallic states shown in 
Fig. [Tj A key characteristic of the deconfined critical point is parametrically large separation 
between the spin correlation length, £, and the confinement scale, £ CO nf, at which hedgehogs 
proliferate. In the doped system, as argued in Ref. [321 the fractionalized excitations are 
already bound into gauge neutral excitations at the scale £. Among these bound states are 
electron-like fermions which carry charge — e and spin 1/2, and so can be directly detected 
in photoemission experiments. We will present a phenomenological effective Hamiltonian 
for these excitations in Section IHIl 



II. U(l) GAUGE THEORY 

The basic ideaP^ of the mapping to the theory of the ACL is to transform to a new set 
of fermions, ipi P with p — ±1, with their spin components p polarized along the direction of 
the local SDW order. We perform this rotation with the spacetime dependent SU(2) matrix 
R av so that™ 

We choose R ap so that spin-fermion coupling is only along <r z , and so 

n«R%a a a ^ p , = a z pp ,=p5 pp , (2.2) 
This relationship is equivalent to 

nl = ^Tr {a a R i (j z R^) (2.3) 

Now, we parameterize 

jt= I " ) (2.1) 




with \ z ia\ 2 — 1- We can verify that Eq. (2.3) yields the usual relation in Eq. (1.2). 
A crucial feature of the resulting Hamiltonian for the ipi P and z ia is that it is invariant 



under a local U(l) gauge transformation. This follows from the invariance of Eqs. (2.1) and 



(1.2) under the transformation 

y. — i. y. p L "<- 



A P -> ^e-**' (2.5) 
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where i?j has an arbitrary dependence on space and time. Note that the ijji P carry opposite 
charges p = ±1 under the U(l) gauge transformation (which is unrelated to the gauge 
invariance associated with the physical electromagnetic force). Associated with this U(l) 
gauge invariance, we will introduce an internal dynamical gauge field in constructing the 
effective theory. 

Ref. |37J argued that describing the transition to the large Fermi surface Fermi liquid 
required the inclusion of additional degrees of freedom so that the theory had a SU(2) 
gauge invariance. We will not consider this extension here. However, we expect that simple 



extensions of the results in Section III apply also to the SU(2) ACL phases found in Ref. [371 



We can now insert Eqs. (2.1) and (1.2) into Eqs. (1.1) and obtain the theory of fluctuating 



Fermi pockets. We will assume that the Zi a are slowly varying, but allow the fermion fields 
ipi P to have an arbitrary dependence on spacetime. The complete Lagrangian is written as 



(2.6) 



The first term is the CP 1 model for the z n : 



C z = - [\(d T - iA T )z a \ 2 + v 2 \(V - iA)z a \ 2 } 
9 



(2.7) 



The fermion hopping term in Eq. (1.1) yields some interesting structure; it can be written 

as 



-E 

Kj 



t%3 



Now, from the derivation of the CP 1 model we know that 



(2.8) 



7* 7 ■ 5» p iJ ^ij 



(2.9) 



and this is easily incorporated into the first two terms in Eq. (2.8), yielding terms which 



are gauge invariant. We therefore incorporate the first two terms in Eq. (2.8) into the 
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gauge-invariant Lagrangian 



p=±l i 



X>U^ + i P A T -ii- AK-lf*^)^ 



(2.10) 



p=±l i<j 



For = 0, £^ describes the band structure in terms of the Fermi pockets. The interactions 
arise from the minimal coupling to the gauge field. Finally, we need to consider the last 
two terms in Eq. ( 2.8 ) . These are the analog of the 'Shraiman-Siggia' couplings.®! Combining 



these terms with the analogous terms arising from the time derivative of the c a , we obtain 
to leading order in the derivative of the z a : 



k,p,q . 



P 



de(k) 
dk 



z±(q - p/2)^(q + p/2)^l(k + q)^+(k - q) + c.c. 



+ ^2(z^d T zn - z il d T z^)ipl_ip i+ + c.c. 



(2.11) 



where e(k) is the single particle dispersion of the large Fermi surface state: 

e(k) = -E^ eik ' (rj_ri) - 



(2.12) 



Note that the terms in C ss mix fermions with different charges. 



The analysis in the following Section III will be based largely on symmetry, and so 
it is useful to recalP^ now how the fields introduced so far transform under symmetry 
operations. These are summarized in Table [T| 



III. EFFECTIVE THEORY OF ELECTRONS AND PHOTONS 

As discussed in Section [I], we want to work in the regime where the photon creates bound 
states between the tp fermions and the z a spinons, but the monopole induced confinement 
has not yet occurred. Thus we are in a fluctuating SDW state with a spin correlation length 
£, but we are interested in phenomena at a scale larger than £. However, the confinement 
of the photons occurs at a scale £ CO nf, and so we will restrict ourselves to the £ < r < £ con f- 

We also note the complementary considerations in the work of Wen and Lee^ they 
considered a spin liquid model with fermionic spinons and bosonic holons (in contrast to 
our bosonic spinons and spinless fermion charge carriers), and described spinon-holon bound 
states in the electron spectral functions. Also, Essler and TsvelilP^S described a model of 
weakly coupled chains, and considered the bound states of the spinons and holons of the 
one-dimensional spin liquid on the chains. In our approach, we do not appeal to these spin 
liquid states, but deal instead with states motivated by the fluctuations in the observed spin 
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TABLE I. Transformations of the lattice fields under square lattice symmetry operations. T x : 
translation by one lattice spacing along the x direction; ' ^0° rotation about a dual lattice 
site on the plaquette center (x — > y, y — > —x); l£ ual : reflection about the dual lattice y axis 
[x —> —x,y — > y); T: time-reversal, defined as a symmetry (similar to parity) of the imaginary 
time path integral. Note that such a T operation is not anti-linear. The transformations of the 
Hermitian conjugates are the conjugates of the above, except for time-reversal of fermiona^H For 
the latter, ip± and are treated as independent Grassman numbers and T : ip± —> —ip±; similarly 
for F a , G a . 

density wave order. 

The regime £ < r < £ con f was treated in Section IV of Ref. [331 See also the subsection 
on holon-spinon binding in the Appendix of Ref. [3U Here we shall provide a more general 
treatment, which should also allow for a computation of spectral functions. 

Let the bound state between the ip + fermions and the z a be F a . The bound state should 
have the full symmetry of the square lattice, and so we can define a local operator which 
creates this bound state centered at the lattice site i. 

As was emphasized in the initial analysis^, there is a second independent bound state, 
and a consequent doubling of the fermion species. This is the bound state between the ifi- 
and z* a1 which we denote by the local operator G ia . In the ordered Neel state, the sublattice 
location of a fermion also fixes its spin. However, when we move to length scales larger than 
£, this is no longer the case because the spin direction of the background Neel state has 
been averaged over. Thus we can view F a and G a as fermions that reside preferentially (but 
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not exclusively) on the two sublattices, and they separately have an additional degeneracy 
associated with carrying S = 1/2. 

More formally, all we will really need are the properties of F a and G a under the square 
lattice symmetry operations: these are summarized in Table |TJ and will form the bases of 
our analysis below. 

The bare electron operator will have a non-zero overlap with both the F a and G a fermions. 



This will be non-local over the scale £. We approximate this connection from Eq. (2.1) as 



^Z(F ta + G lct ) (3.1) 

where Z is some quasiparticle renormalization factor depending upon the fermion-spinon 
bound state wavefunction. In general, Z should be non-local over a scale £, but have limited 
ourselves for simplicity to a momentum independent wavefunction renormalization. Note 
that Eq. (3.1 ) and the symmetry transformations in Table [T] ensure that c a is invariant under 



all operations of the square lattice symmetry. The possible non-local terms in Eq. (3.1) can 
be deduced by the requirements of symmetry. 

We now need an effective Hamiltonian for the F ia and Gi a . Formally, any Hamiltonian 
which is invariant under the symmetry transformations of Table [I] is acceptable; however, we 
use simple physical requirements to restrict the large class of possibilities. For the diagonal 
terms which do not mix the F and G, we assume (for simplicity) that they just inherit the 



terms for ip + and "0- in m ^1- (2-10). The mixing between the F and G is provided by the 
Shraiman-Siggia term C ss . Physically this can be understood as the mixing corresponds to 
hopping between two sublattices (as F and G reside preferentially on different sublattices), 
and Shraiman-Siggia term describes such hopping and associated spin-flipping process®!. 
These terms are more simply considered in their real space form which are the last two 
terms in Eq. ( |2.8[ ). Combining these terms, we have the effective Hamiltonian 

H cS = - ^ Uj (F^Fja + G\ a Gja) - A ^2(-l) lx+l y (^F- a F ia - G\ a G ia j 

ij i 

-}Z^{ f L g ^ + g L f j«) (3-2) 

where the second line comes from C ss . Here tij and tij are renormalized in some unknown 
manner from the bare hopping matrix elements in C s f. We can now verify that all the terms 
above are invariant under the transformations in Table HI 

The above Hamiltonian of bound state can be diagonalized in momentum space. First, 
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rewrite equation ( |3.2[ ) in momentum space 

H cS = J>(k) - M ) (i^F*. + GL.Gte) + J>(k) + G^fJ 

~ (^ka-^k+K,a - G^Gk+K^j (3.3) 

where we parameterize e(k) and e(k) as 

e(k) = — 2t(cos k x + cos k y ) + At' cos fc x cos k y — 2t"(cos 2fc x + cos 2k y ) (3.4) 
e(k) = —t — 2t(cos k x + cos k y ) + 4? cos k x cos fc^ — 2t"(cos 2k x + cos 2^) (3.5) 

Here t, if and t" are nearest neighbor, next nearest neighbor and next-next-nearest neighbor 
hopping tij respectively, t, t' and t" are the hopping elements of Uj. Here we only included 
terms up to third nearest neighbor hopping, which is capable to capture the shape of the 
Fermi surface, but higher order terms can be included in a similar fashion. i is the matrix 
element of the on-site mixing term F}Gi + G\Fi, which is also allowed by symmetry. 



To diagonalize Hamiltonian (3.3), we change basis to 



C*k« — y|^ kQ ^ kcf )' — ^^(Fk+K,a — Ck+K,a) (3.6) 

Note that D^ a has momentum k according to the transformation under translation symmetry 
listed in Table [TJ In the new basis, the Hamiltonian becomes 

H cS = [(e( k ) + - v)Ct a C ka + (e(k + K) - e(k + K) - fi)Dt a D ka 

k 

-A^L^ka + ^LCka)] (3.7) 

The spectrum of electron operator c can be obtained by diagonalizing the Hamiltonian 



of C and D fermions, as c is related to C operator according to equation (3.1) c^, 



(Z/y/2)C ka . In the Hamiltonian (3.7), C and D fermions have dispersions e(k) + e(k) 
and e(k + K) — e(k + K) respectively, and they are mixed through the A term. With the 
mixing, gaps will open where the Fermi surfaces of C and D fermions intersect and the 
large Fermi surfaces become Fermi pockets. In the case of e(k) = 0, the D Fermi surface 
is the same as the C Fermi surface shifted by so the pockets are symmetric under 

reflection with respect to the magnetic Brillouin zone boundary, and therefore centered at 
(7r/2,7r/2). However, with a non-vanishing e(k), the dispersions of C and D are different, 
so the pockets are no longer symmetric about the magnetic Brillouin zone boundary, and 
are not necessarily centered at (tc/2,it/2). 

To show the qualitative effects of e on the shape of Fermi pockets, we draw the pockets and 
electron spectrum functions of some representative choice of e in Figs. |2](5j The dispersion 



e(k) is chosen with some phenomenological parameters if = 0.15£ and t" = — 0.5£', and the 
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FIG. 2. Plot of Fermi pockets of the bound state with e(k) = 0.8i. Subplot (a) shows the 



Fermi surface of Ck Q and -Dka fermions as if there is no mixing term A in the Hamiltonian (3.7). 
The red line is Ck« and the blue one is D^ a . Subplot (c) shows pockets like Fermi surface of the 



quasiparticles described by the Hamiltonian (3.7) with A = 0.3i. Subplot (b) shows the same Fermi 



surface with the color representing the weight of electron operator in the quasiparticle excitation. 
Subplot (d) shows the same information as in (b), but by plotting the electron spectral weight at 
lo = as a function of momentum; in this plot we manually put in a finite life-time of electron 
r ~ 0.2t just for visualization purpose. The dashed line in (a) and (c) is the boundary of magnetic 
Brillouin zone in the ordered state. 



SDW gap A is 0.3t in Figs. |2|4]and 0.5t in Fig. [H] In Fig. [2J a negative t shifts the hole 
pockets outwards, and makes the shape of the hole pockets asymmetric. The inner side 
becomes more curved and the outer side more flat. In Fig. [3J a negative t does not shift 
the position of the hole pockets significantly, but also makes the shape of the hole pockets 
asymmetric in a similar way as in Fig. [2j Combining the effect of these two parameters, 
we can move the hole pockets inward and make their inner side more curved than the outer 
side with a position to and a negative t, as shown in Fig. |4j Along the Fermi pockets, 
the fermionic quasiparticles are a mixture of Ck Q and Dk a fermions, while experiments only 
probe electron spectrum weight. The weight of electron operator in the quasiparticle is 
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-7T -7r/2 T/2 7T 



FIG. 3. As in Fig. 15] but with e(k) = 0.5t(cos k x + cos k y ) and A = 0.3t. 



calculated through diagonalizing equation (3.7) and are plotted in Fig. [2j|4j In Fig. |5| we 
show a plot with larger A (A = 0.5i) so that the anti-nodal electron pocket is completely 
gapped. In all the cases, the inner half of the hole pockets have higher electron quasiparticle 
weight, since the inner part is primarily made of C^q fermion and the outer part is primarily 
made of Dy_ a fermion, and electron operator is proportional to Ck Q . 

This bound state theory can be compared to the Yang, Rice, and Zhang (YRZ) model 
of hole pockets^^. In their theory a phenomenological Green's function for the electron 
in the underdoped state is proposed, based upon "spin liquid" physics, and the spectral 
function derived from the Green's function has hole pockets inside the diamond Brillouin 



zone. From our Hamiltonian (3.7), the Green's function of the electron is 



G c (k,co) 



u - e(k) - e(k) + n + X 2 /(u - e(k + K) + e(k + K) + /i) 



(3.8) 
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FIG. 4. As in Fig. [2] but with e(k) = -0.3i + 0.5i (cos fe x + cos and A = 0.3*. 

Compare this to the Green's function in the YRZ modeP^E 

ft 



G(k,u) 



(3.9) 



The two results are quite similar. In the present form, instead of the ci-wave gap our A 
does not have a momentum dependence. However, this distinction is an artifact of the simple 
choices made in our form of H e g: we can clearly include more non-local mixing between the 
F and G fermions. 

The more important distinction between our model and YRZ lies in the physical input 
in the fermion spectrum. The YRZ model relies on 'pairing correlations' implicit in some 
underlying spin liquid state. In contrast, we do not assume any pairing, but fluctuating 
local antiferromagnetic order. We will describe the influence of pairing on our spectral 
functions in a future paper. We note earlier arguments^ that the pairing amplitudes are 
especially strong on the electron pockets, and this may be a contributing factor to removing 
the electron-like Fermi surfaces in Figs. [3] and |4| 

Recent ARPES experiments^ reveal that there are hole pockets in the underdoped regime, 
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FIG. 5. As in Fig. [2] but with e(k) = -0.3* + 0.5t (cos fe a + cos and A = 0.5i. 



and the center of the pockets are inside the first magnetic Brillouin zone. In general the 
shape of hole pockets seen in the experiments can be fit to our model. As shown in Fig. [2|j5j 
the weight of electron operator on the outer part of the pockets is tiny and may be hard to 
see in experiments. 

We close this section by remarking on the status of Luttinger's theorem in our theory 
of the electron spectral function. These issues were discussed in Refs. [3U and H3j where 
we argued that the theorem applied to the sum of the F a , G a and the if)± Fermi surfaces. 
Here we will focus on the electron-like F a , G a , and will drop the ip± contributions to the 
present discussion; the latter will amount to shift in the effective doping level x. Under this 
assumption the total number of F a and G a per site is 2 — x, and so is the number of C a 



and D a , since the canonical transformation in equation (3.6) preserves particle number. On 



the other hand, our theory was applied to the doped holes or electrons in the background 
of a fluctuating antiferromagnet, and each such charge carrier must occupy one state within 
the Fermi surface. Counting hole (electron)-like Fermi surfaces as negative (positive), then 
for a doped antiferromagnet with hole density x, the Fermi surfaces in Figs. |2]j5] should 
enclose a total area of — (27r) 2 (x/2) (the last factor of 2 is from spin degeneracy). In our 
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present effective Hamiltonian model, we have found it more efficient to treat not just the 
fermionic excitations near the electron and hole pockets, but across the entire Brillouin 
zone. We found that this method was very convenient in treating the contraints imposed by 
symmetry, without prejudicing the final locations of the Fermi surfaces. Such an extension 
should not be accompanied by any fundamental change in the many body quantum state, 
and hence cannot modify the statement of Luttinger's theorem. Because the Luttinger 
constraint only controls the electron density modulo 2 per unit cell, we therefore conclude 
that our lattice effective model obeys 

F} a F ia ) + (G\ a G ia ) = (cj a C ia ) + (plD ia ) = 2 - x (3.10) 

Note that this differs from the value in the conventional Fermi liquid phase, in which case the 
total electron density is 1 — x. This difference is acceptable here because we are discussing 
a phase with topological order, which has an emergent U(l) gauge excitation®'. 

In the discussion of YRZ modeP^^ a different form of Luttinger's theorem is used: there 



the total number of particles given by Green's function (3.9) is given by the total area 
enclosed by the contours where the Green's function changes sign. This includes not only 
the Fermi pockets on which the Green's function diverges, but also the contour where the 



Green's function vanishes. For our Green's function of C bound states in equation (3.8) 
this extra Luttinger surface of zeros is exactly the original Fermi surface of D bound states 
without the mixing between C and D. Therefore the total area enclosed by both Fermi 
pockets and surface of zeros equal to total number of C and D states minus the number of 
D states, and the result is the total number of C states. So this form of Luttinger's theorem 
still holds true in our model, although it is not relevant to the physical doping, which is 
related to the area of Fermi pockets only. 



A. Coupling to photons 

Apart from the electron-like bound states described in the previous section, the other 
low-lying excitations in the ACL are the gapless U(l) photons. Here we discuss how these 
two low-lying sectors of the theory couple to each other. 

Our first task is search for terms coupling the F a , G a to the photons, while being invariant 
under all the transformations in Table [Tj It is useful to do this on the lattice, as many of the 
operations involve details of the lattice symmetry. The F a , G a are gauge- invariant, and so 
will couple only to the field strengths: on the lattice, we define these as 

B = A x A y — A y A x , £ x = A. X A T — A. T A X , £ y = A y A T — A T A y ; (3.11) 

thus B resides at the center of each square lattice plaquette, while £ x>y reside on the links. 
After some searching, we found a single term which is linear in the field strengths which 



16 



fulfills the needed criteria: 



= -z 7 fdrJ2 % (FlKG a - GlA x F a ) + £ y [FlA y G a - G* A y F a )] (3.12) 

i 

This coupling contributes to the self-energy correction of the bound states. At the lowest 
order, the self-energy correction is represented by the Feynman diagram in Fig. [6j 




FIG. 6. Feynman diagram for the self-energy correction. The solid line with arrow represents F 
or G fermion propagator and the wiggly line represents electric field propagator. The interaction 



vertex is given by equation (3.12). This diagram is evaluated in equation (3.18). 



The action of U(l) gauge field has contributions from both the ip± fermions and z a field. 
Assuming that there exist Fermi surfaces of the ip± fermions, integrating them out screens 
the fluctuation of the A T component, and gives the following terms for the action of the 
transverse components of A 

Here k F is the Fermi momentum of the ip± Fermi surface. 

We consider two possibilities of the z a fields. At the magnetic ordering critical point, 
integrating out z a field gives rise to the following action 

s Az = \j j0M-p, -<*) (f - ^f) V^^Mp,^) (3-14) 

On the other hand, if the z a field is gapped as in the ACL, integrating it out gives rise to 
terms proportional to u 2 + cj,p 2 , which are higher order than those in Eq (3.13). 
In summary, the U(l) gauge field action can be written as 

sAz = U ^M-p,-^MpMMp^) (3-15) 

where the polarizition function II is given by 

II(p, uj) = kp^- — - + a/w 2 + c 2 p 2 + ■ ■ ■ at the quantum critical point (3.16) 

II(p, u) = kt— + P 2 in the ACL (3.17) 

P 
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for imaginary frequencies u. 

Using the action in equation (3.13) and (3.14) we can evaluate the Feynman diagram in 
Fig. [61 



E(p,w) = 7 2 



d 2 qdtt 



1 



n 2 



(2tt) 3 i(n + u)-^ q n(\p-q\,Q) 



{Pi + qi)(pj + Qj) 



5 ij 



iPi- Qi)(Pj- Qj) 

- qf 



(3.18) 

Here we consider a simple one-band model for the bound state with a quadratic dispersion 
£ a = q 2 / (2m) — fi. The polarization factor at the end can be simplified as 



(Pi + Qi)(Pj + Qj) 



8 13 - 



iPi ~ Qi)(Pj - Qj) 

(p - q) 2 



Ap z q 2 (l - cos 2 9) 
p 2 + q 2 — 2pq cos 6 



(3.19) 



where 9 is the angle between p and q. 

We are interested in the life-time of quasiparticles near the Fermi surface, so we consider 
the imaginary part of self-energy at the quasiparticle pole (p,cu = £ p ). Without losing 
generality, we consider the case of u > 0. At zero temperature, energy of the intermediate 
fermion state must be greater than zero but smaller than oj, so its momentum must be inside 
the shell of kp < q < p. The imaginary part of self-energy is 



ImS rct (p,w) 



r 



(2tt) 2 
Im 



2?r 



qdq / d9(u — £ q 



Ap 2 q 2 (l - cos 2 9) 



(3.20) 



Il ret (|p — q\, uj — £ q ) p 2 + q 2 — 2pq cos 9 



Near the Fermi surface we can linearize the dispersion relation £ q = Vp(q — Pf), and 
uj = £ p = vp{p — Pf). We change the integrated variable from q to k = p — q, < k < lo/vf- 
In the limit oi u — > 0, oi p — pp — > 0, the integral can be simplified 



ImS ret (p,w) 



f]4 I v f 
7T 2 m 2 



2- 



dk / d9k 2 lm 



1 + cos 9 



U^(^k 2 + 2p 2 F {l -cos 6), v F k) k 2 /{l - cos9)+p 2 F 

(3.21) 



In both of the cases considered in Eqs. (3.16) and (3.17), the gauge field spectrum function 

1 



Im 



W et (^k 2 + 2p 2 F (l-cos9), v F k) 



AT 1 , A;^0 



when 9^0, and becomes less singular at 9 = 0. As a result the integral in Eq. (3.20) has 
the following behavior 

Im£ ret (p,w) ~ oo 2 , oo^O. 

Thus the bound state fermion has a Fermi liquid-like damping. 

This conclusion is different from Essler and Tsvelik's work on a model of weakly coupled 
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chaindHEH. Their pocket states arose from binding between the one-dimensional holons 
and spinons, and the coupling to the one-dimensional spin fluctuations led to a self-energy 
proportional to a;ln(a;). Our model is genuinely two-dimensional, and has an emergent 
U(l) gauge field who presence is also the key to the violation of the Luttinger theorem in 
Eq. (3.10). They explain the violation of the Luttinger theorem by zeros of the Green's 
function 1 54 1 55 1 58 1 ; the connection of these zeros to our work was discussed below Eq. (3.10). 



IV. ELECTRON SPECTRAL FUNCTIONS IN THE RENORMALIZED CLASSI- 
CAL REGIME 



We will now consider the regime at small non-zero T above the antiferromagnetically 
ordered state present at g < g c , the renormalized classical (RC) regime. Here, we expect 



the 0(3) vector formulation of C s f in Eq. (1.1) to be equivalent to the spinor/U(l) gauge 
theory formulation of C ac i in Eq. (2.6). This is explored at T = in Appendix [Al where we 



show that the small g expansions of the two theories match with each other. Thus, while 
there are differences in the nature of the non-magnetic phases at g > g c in the two cases (as 
reviewed in Section [T]), the g < g c phase and its low T RC regimes are the same. 

The RC regime was studied using the 0(3) formulation by Vilk and Tremblay 49 , and 
using a U(l) gauge theory similar to C ac i by Borejsza and Dupuis^l. Here we shall also use 
C ac i, and expand upon these earlier results. 

We compute the electron spectral function as a convolution of the z a propagator with 
the free if) propagator containing pocket Fermi surfaces obtained from C^. 

For the z a propogator we use the simple damped form motivated by studies in the l/N 
expansion.^ 

G z {k,u) = — — — r^- — 2 (4.1) 

where the dimensionless constant 7 determines the damping constant, and the "mass" m is 
determined by the "large N" equation 



m = IT In 



(4.2) 



where 

1 1 

V S 9c 

is the energy scale which determines the deviation from the quantum critical point. The 



Q 



(4.3) 



damping in Eq. (4.1) is a simple interpolation form which is constant with the expected 



behavior at the quantum critical point, and with the RC behavior of Eq. (7.28) in RefP} 
Actually, we have neglected the power-law pre-factor of the exponential in (7.28). Determin- 



ing the damping from numerical studies motivated by the e expansion^ 1 we expect 7 fa 1. 
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In Fig. [7] and [8] we plot electron spectrum function obtained as a convolution of the z a 
and fermion spectrum functions at finite temperature. Here fermion dispersion relation is 



the same as the one we used in the plots in section III In contrast to Fig. [2}f5l the pockets 



are symmetric with respect to the magnetic Brillouin zone boundary, and hole pockets are 
centered at (n/2,7r/2). Again, the inner half of the pockets has higher spectrum weight, as 
the quasiparticle on outer half of the pockets is primarily made of fermion quasiparticles at 
fk+(Tr,n)- The broadening of the spectrum comes from two factors: first, the z boson has a 
finite life-time 7m; second, the convolution is done at finite temperature, so the energy of 
fermion and boson excitation does not have to be exactly zero, but can vary by energies 
of order T. The second factor dominantes in the RC regime because m < T, and 7 is of 
order 1. So the electron spectrum is basically the pocket fermion dispersion broadened by 
a linewidth of order T: see Fig. [7] and Fig. [8j Earlier workP^ had also obtained a linewidth 
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FIG. 7. Electron spectrum function at zero frequency as a convolution of z a and the ip± fermions. 
In the plot we chose the following parameters: g = O.lt, T = 0.2t and the SDW gap A = 0.2t. The 
white dashed line shows the boundary of the diamond magnetic Brillouin zone of the commensurate 
(tt, 7r) AFM order. 



of order T in the RC regime by very different methods. 

As we noted above, the present Fermi surface locations are symmetric with respect to the 
magnetic Brillouin zone boundary, unlike those in Figs. [2]{5j However, as shown in Ref. 



the Shraiman-Siggia term in Eq. (2.11 lifts this symmetry. We can compute perturbative 



corrections to leading order in A, and these will contribute spectral weight of width T which 
is asymmetric about the magnetic Brillouin zone boundary. 
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FIG. 8. Electron spectrum function at zero frequency as a convolution of z a and the ip± fermions. 
In the plot we chose the following parameters: g = O.lt, T = 0.2t and the SDW gap A = 0.5t. The 
white dashed line shows the boundary of the diamond magnetic Brillouin zone of the commensurate 
(-7T, 7r) AFM order. 

V. CONCLUSIONS 

We conclude by reviewing the different routed by which the spin fermion model C s f 



in Eq. (1.1) can lose antiferromagnetic order, and their implications for the photoemission 
spectrum at finite temperatures. 

The first, more conventional, route is that there is a single direct transition at g = g c to a 
Fermi liquid with a large Fermi surface. Then in the T-g plane, we have the conventional 
quantum disordered (QD), quantum critical (QC), and renormalized classical (RC) regions. 
In QD region at small T for g > g c , the electron spectrum will show quasiparticle peaks 
with a Fermi liquid linewidth ~ T 2 . In the QC region near g = g c , the spectrum will 
again show weight along the large Fermi surface, but with large anomalous linewidths near 
the 'hot spots': these are points along the large Fermi surface connected by the ordering 
wavevector K. Finally, in the RC region at small T for g < g c , we have the behavior 



described in Section IV the spectrum has 'small' Fermi pockets which are centered at the 
antiferromagnetic Brillouin zone boundary, and the quasiparticle peaks have a width ~ T; 
examples of such spectra were shown in Figs. [FJ |8} 

This paper mainly considered a more exotic route^Myj towards loss of antiferromagnetic 
order in the spin-fermion model. This route is possible if topological "hedgehog" tunnelling 
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events are suppressed as the transition, as they are at magnetic-disordering transitions in 
the insulator (see Fig. [T|. Then the transition at g = g c is to a non- magnetic non- Fermi 
liquid ACL, with the spinless fermions ip± and complex bosonic spinons z a as elementary 
excitations, interacting (in the simplest case) via an emergent U(l) gauge force. Again, 
around this critical point at g — g c , we can define the corresponding RC, QC, and QD 
regions. The electron spectrum in the RC region is just as in the first case above, as 



discussed in Section IV and Appendix [A] The novel spectrum in the QD region was the main 
focus of Section III One contribution to the spectral weight comes from the convolution of 
the deconfined spinons and the ip±. This leads to incoherent spectral weight in 'arc'-like 
regions which were described earlier in Ref. [331 However, it has been arguecP^' that the 
spinons and ip± form electron-like bound states, and these were described in more detail 



in Section III We found that the bound states lead to pocket Fermi surfaces, as shown in 
Figs [2](5] An important feature of these spectra are that the pockets are not centered at 
the point (7r/2,7r/2) on the antiferromagnetic Brillouin zone boundary. In fact, this zone 
boundary plays no special role, and there are no symmetry relations on the quasiparticle 
dispersions across it. We note that recent photoemission observations^ show features related 
to this QD region. We did not address the QC region here, but it should mainly have the 
incoherent arc spectra, as described in Ref. 1331 

57H591 



We also noted a similarity between our ACL QD results, and the YRZ mode l " ' in 



Section |III| The assumptions of YRZ are very different from ours, as they depart from a 
paired fermionic spinons in a spin-liquid state. The Fermi surfaces in the YRZ model to 
not obey the conventional Luttinger theorem, as is the case in our model. However, such 
a violation must be accompanied by gauge forces reflecting the topological order in such a 
stated these gauge fields do not appear in their formulation. Also, we dealt mainly with the 
influence of local antiferromagnetic order on the electron spectrum. We have not included 
pairing effects in our computations yet, or the transition to superconductivity; these are 
issues we hope to address in forthcoming work. 
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Appendix A: Small g expansion 



This Appendix will compare the properties of the Lagrangian C s f defined in Eq. (1.1), 



with the 'fractionalized' Lagrangian C ac i in Eq. (2.6). We will work at T = in the limit of 



small g, where both models have long range SDW order, and are expected to be essentially 
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identical. We will compare the two models here by computing the on-shell matrix element 
for the decay of a spin-wave in the SDW state into a fermionic particle hole pair. 

We will work in the "diamond" Brillouin zone associated with antiferromagnetic ordering 
with wavevector K = (ir, n). In this zone we define 

ci Q (k) = c a (k) , £i(k) = e(k) 
c 2a (k) =c a (k + K) , e 2 (k)=e(k + K) (Al) 

All expressions below are implicitly in this diamond Brillouin zone. 

The analysis of the spin-wave decay appears for the two models in the following subsec- 
tions. 



1. Spin-fermion model 

We perform the small g expansion for the order parameter by the following parameteri- 
zation in terms of the spin-wave field 0: 



(t^L^g-gW*, ^fV2g-9 2 \<P\ 2 , 1 - 9\<P?) (A2) 



Then the Lagrangian for n a is 

C n = Id^l 2 + f [(0*<9 M 0) 2 + (0<V>*) 2 ] + y (2 |^ |2) ( W + <pd^f (A3) 

and can be analyzed as usual in an expansion in g. For the fermion sector, we diagonalize 
the C c + C\ at g = by introducing fermion operators 7i,2 P , and (we replace the electron 
index a by p = ±1) 

cip(k) = M k 7ip(k) - jWk72 P (k) 

c 2p (k) = 2w k 7ip(k) + Mk72 P (k) (A4) 

where «k, t>k are real and obey w 2 , +v%. = 1. We choose = cos(#k/2), t>k = sin(0k/2), and 
then 

£l (k)-e 2 (k) . ~2A 
cost^k = — sm# k = — . A5 

V /(£ 1 (k)-£ 2 (k))2+4A 2 V /(£ 1 (k)-£ 2 (k))2 + 4A 2 

The fermion Lagriangian at g = is 

£7 = E 7k & + ^i( k )) Tip + E ^1 & + ^00) 72 P (A6) 

k,p k,p 
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where 



E 1>2 (k) = - L(k) + e 2 (k) ± v / (^i(k)-e 2 (k)) 2 + 4A 2 



(A7) 



Finally, the non-linear couplings between the spin waves and the fermions is given by 



k,q 



J q 



(wk+q^k - n k v k+q )(7 1 t _(k + q)7i+(k) - 7j_(k + q)7 2 +(k)) 
+ Ok+qWk + u k u k+q )(7l_(k + q)72+(k) + 7L( k + q)7i+(k)) 



H.c. 



+gX 0^| 2 ] q K+q^k + n k i; k+q )(7 1 t p (k + q)7i„(k) - 7 2p (k + q)72 P (k)) 

k,q p 

+ p(li k+q M k - Wk'^k+q )( 7 l p (k + q) 72p (k) + 72 t p (k + q) 7 i P (k)) (A8) 
Now we can obtain the self-energy of the <fi spin-wave excitation to order g: 



(^k+q^k - ^k^k+q 



+ (u k+q u k + f k W k+q 



i2/ /(£ 1 (k))-/(£; 1 (k + q)) /(£ 2 (k))-/(£ 2 (k + q)) 



E 1 (k + q) - E 1 (k) -iu n + E 2 (k + q) - E 2 (k) 
f{E 2 {k))-f{E x {k + ^)) /(E 1 (k))-/(E 2 (k + q)) 



+ 



■iu n + E 1 (k + q)- E 2 (k) -iu n + E 2 (k + q) - E^k) 
+ 2g\J2^ k v k (f(E 1 (k)) - f(E 2 (k))) 



(A9) 



As expected, we have E^(0, 0) = 0. We can also estimate the on-shell decay rate at T = 0: 
ImE^fvg, q) ~ gg 2 , as long as v < Vf- 

We can also use to compute the fermion spectral density at order g. For the case of 
a single particle in the insulator, it seems to me that the result has the same form as (2.47) 
in Ref . 



2. U(l) gauge theory 

For small g, we parametrize the z a field of the CP 1 model as 



-id 



(AlO) 
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Then the Lagrangian for z a is 



\du 



+ 



4(2 -gW) 



(All) 



If we integrate over An in the gauge i? = 0, we reproduce the action in Eq. (A3). 



The integral over A^ also produces a coupling between 0*6^0 — 0<9 M 0* and the U(l) 
current over the ipi^p fermions. Note that this is a quartic coupling, a bilinear in <fi coupling 
to a bilinear in ipi,2p- 



In our discussion in Section A 1 , the dominant coupling between the fermions and the 
spin-waves was the term in L$ a which was linear in <fi. Such terms arise here exclusively 
from C ss . We evaluated the on-shell matrix elements from such terms for the decay of a 
spin wave into a particle hole pair of the ji p (and also to the j2p)- From the first term in 



C ss we have the matrix element 



q 



dei(k) 2 <^2(k) 

~dk~ ~ Vk ~dk~ 



(A12) 



For the time-derivative term in C ss we replace the frequency by its on-shell value q ■ 
dE\ (k) / <9k to obtain the matrix element 



(A13) 



The sum of (A12) and (A13) is equal to the matrix element obtained in Section Al which 
is 

(A14) 
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